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Small-time asymptotics for Gaussian self-similar stochastic 

volatility models 

Archil Gulisashvili,* Frederi Viens,l and Xin Zhang,1 


Abstract 

We consider the class of self-similar Gaussian stochastic volatility models, and compute the 
small-time (near-maturity) asymptotics for the corresponding asset price density, the call and 
put pricing functions, and the implied volatilities. Unlike the well-known model-free behavior 
for extreme-strike asymptotics, small-time behaviors of the above depend heavily on the model, 
and require a control of the asset price density which is uniform with respect to the asset price 
variable, in order to translate into results for call prices and implied volatilities. Away from 
the money, we express the asymptotics explicitly using the volatility process’ self-similarity pa¬ 
rameter H, its first Karhunen-Loeve eigenvalue at time 1, and the latter’s multiplicity. Several 
model-free estimators for H result. At the money, a separate study is required: the asymptotics 
for small time depend instead on the integrated variance’s moments of orders | and |, and the 
estimator for H sees an affine adjustment, while remaining model-free. 
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1 Introduction 

In this paper, we present a study of the small-time (near-maturity) asymptotics for the asset price 
density S, the call and put prices, and the implied volatilities, for the class of continuous-time Black- 
Scholes-Merton-type models with Brownian noise and independent Gaussian self-similar volatility. 
The techniques borrow from a framework established in our prior work |44j for general Gaussian 
volatility models; they use a tailored application of Laplace’s method requiring a delicate analysis of 
uniformity with respect to strike prices K away from the money {K sq), and apply a general result 
from [35] to translate asymptotics from call prices to implied volaltilities. Model-free estimators of 
the self-similarity parameter H result. Away from the money, all asymptotic constants and powers 
are expressed explicitly in terms of H and of the coefficients in the Karhunen-Loeve expansion of 
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the volatility. At the money {K = sq), a separate study is required. This introduction contains 
extensive details of general context of the small-time asymptotic problems mentioned above, our 
motivations, and a precise summary of all our results. 

1.1 General background 

It has been known for decades that the Bachelier-Black-Merton-Scholes framework, while extraor¬ 
dinarily fertile for explaining various basic features of financial markets and for helping define 
fundamental notions, including volatility as the relative scale of noise intensity, suffers from certain 
deficiencies, particularly the fact that volatility is not constant empirically. When coupled with 
the fact that non-random volatility, which implies normally distributed log returns, has difficulties 
in explaining certain extreme events because of excessively light tails, one quickly arrives at the 
vast class of stochastic volatility models, i.e. those continuous-time models where the relative noise 
intensity of returns is itself a stochastic process which is at least partially driven by exogenous 
noise. A large number of articles and monographs on stochastic volatility (SV) can be consulted for 
empirical and economic justification of these models; we cite the classical text [32]. Of particular 
interest is SV models’ ability to reproduce some desirable market features of option prices, such as 
“smiles” and other non-flat shapes of the implied volatility (IV), i.e. the volatility which would be 
required of a constant-volatility model to explain a given call option price. 

One of the first mathematical treatments explaining empirically observed IV shapes was by 
Renault and Touzi in [SS]. Recent studies have looked in detail at the question of IV asymptotics, 
that is to say the behavior of IV as important parameters such as strike price K and maturity T 
tend to extreme values. Of note is the groundbreaking paper |49j of Lee, in which the large-strike 
(the small-strike) behavior of IV is described in terms of the largest (the smallest) non-exploding 
moment of the stock price. Gaussian volatility models belong to the class of models with moment 
explosions. For more details and other references on IV shapes and extreme-strike asymptotics 
of IV, we refer to the introduction section in our prior work [44] . where we examine the class of 
uncorrelated Gaussian volatility models in its broadest possible sense. 

1.2 Specific motivations and modeling choices 

Small-time asymptotic behavior of densities, option pricing functions, and implied volatilities has 
been a popular topic of study. There are various model-independent results (see, e.g., 01331171 IMj), 
explaining how the asymptotics of the IV depend on those of option pricing functions. There are also 
papers discussing small-time asymptotics of the functions mentioned above in the case of stochastic 
volatility or local-stochastic volatility models (see IlllTlEaETlEHlIMl sa ETj), and for special 
models (see [3J [HI [33 [HI jSH |S3| (models with jumps), [331 [33 [Ml |3D| (Heston model), [13 [30] 
(Stein-Stein model), [13I13I111IS7] (SABR model)). 

The present paper follows up on our prior study in [H] by attempting to elucidate the small¬ 
time behavior of IV for a subclass of Gaussian volatility models, consisting of models with self¬ 
similar volatility processes. It turns out that establishing small-time asymptotics in a general 
Gaussian context is significantly more demanding than determining large-strike behavior. This can 
be understood as a manifestation of the fact that there is no model-free analogue of Lee’s moment 
formulas in the small or large time regimes. In this paper, we illustrate the challenge by specializing 
to the case of self-similar volatilities; we will see that the type of small-time behavior for both call 
price and IV is quite sensitive to the self-similarity parameter H. This is good news if one is to 
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leverage these results to help determine H, as we will see. 

Indeed, our study also allows us to investigate the question of long-memory SV calibration, since 
long-range dependence and self-similarity are proxies for each other in many known models, via their 
common Hurst parameter H. Based on a Gaussian long-memory model for log-volatility pioneered 
by Comte and Renault in E], the work in uni used an ad-hoc calibration method based on option 
prices to determine H so as to best explain market prices. Fractional volatility models also appear 
inHHaiMUSZlEHlIlIlESIMliniEaEI]- In the current paper, we show that calibration of H 
near maturity can be given a stronger mathematical foundation under self-similarity assumptions 
for the volatility process. The parameter H can also be a proxy for local regularity measurements, 
in the sense of their paths’ Holder continuity parameter. Some recent papers and presentations, yet 
unpublished at the time of writing this article, appear to show that volatility is rough, in the sense 
that the log-volatility process is fractional and it is not Holder continuous for 1/2 — e < H < 1/2, 
where e is a positive number (see [36l [37l |38] ). On the other hand, [10] and many studies before it 
(see references therein) indicate that i? > 1/2 in terms of memory length. This is a demonstration 
that the use of H to measure self-similarity and long memory and path regularity/roughness, such as 
in the case of fractional Brownian motion (fBm), might be a misspecification in volatility modeling. 
The authors of [35] indicate that classical long-memory tests detect this property in their Gaussian 
rough volatility model, which is a geometric fBm or a geometric OU process with shorter memory 
{H < 1/2). The studies in [10] show on the other hand that no consistent memory estimation results 
in practice from any classical method when used on the non-self-similar stationary long-memory 
model of m- Our current work could help in elucidating the differences between these points of 
view; we do not comment on them further herein. An interesting discussion of long memory vs 
short memory problem can be found in Section 1.2 of [118] . In any case, the numerics which we 
include in this paper and will discuss at the end of this introduction show that our model class 
allows for a very sharp calibration tool. 

Before providing a summary of our results, we discuss some classical Gaussian self-similar mod¬ 
els. General details about this class are given in Section [^ These are the Gaussian processes 
X on [0,r] such that for some H G (0,1) and for any a > 0, the two processes t i —Xat and 
t I —Xt have the same distribution (law). The best known among them is the fractional Brow¬ 
nian motion (fBm) , the centered Gaussian process whose law is defined by (0) = 0 and 
E (B/^ — ~ s|^^. It is the only (continuous) self-similar centered Gaussian process with 

stationary increments. Many texts can be consulted on B^, including, e.g., [331133113T] . Among 
the many other centered Gaussian self-similar models, which are all necessarily non-stationary. 


H,RL 


= /o (i - s) 


H-l/2 


dW{s) 


the easiest to construct is the Riemann-Liouville fBm, defined as B, 
where W is a standard Wiener process (see for instance [50]). This process, which is B-self-similar, 
has properties close to those of fBm, and can be more amenable to calculations. The so-called 
Bifractional Brownian motion depends on two similarity parameters H and K, has a more complex 
representation, as the sum of an fBm with parameter HK, and a process with paths which is 
not adapted to a Brownian filtration: see [48], see also [3] and the references therein. This process, 
which is BAT-self-similar, can model the effect of smoothly acquired exogenous information, and 
is an extension of the so-called sub-fractional Brownian motion (see 0). Self-similar Gaussian 
processes can also be obtained as the solutions of stochastic partial differential equations: a class 
which includes solutions to fractional colored stochastic heat equations is studied in [^ , which has 
the interesting property that its discrete quadratic variation has fluctuations which become non- 
Gaussian at a threshold of self-similarity which is lower than for fBm, and can be adjusted to be 
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as low as desired. This can be helpful to model volatilities whose local behavior has heavier-tailed 
fluctuations than what standard fBm can allow, regardless of the volatility’s self-similarity. It also 
allows the modeler to choose regularity and self-similarity properties independently of each other, 
which offers more flexibility than the models considered in [nmniisH]. More examples of Gaussian 
self-similar process can be found in Interestingly, many of the Gaussian self-similar models 


share the same path regularity properties as fBm, because it can be shown that there are positive 


finite constants c,C for which c\t — s 


,2H 


< E 


\Xt — Xa\‘^ < C\t — s'^^, where the symbol H 


stands for the self-similarity parameter of the model under consideration. 

Finally, it bears noting that self-similarity implies that Xq = 0 and that Var [Xt] is proportional 
to . This is a strong assumption on X. An uncertainty level on volatility which increases with 
time is a reasonable conservative forecasting assumption. That the volatilty starts at 0 is more 
restrictive, since, in our IV context, it corresponds to saying that the underlying risky asset’s 
movements tends towards certainty near the derivative’s maturity. Such a behavior is characteristic 
of specihc risky asset classes, such as hxed-income securities, e.g. treasury bonds, and the dividend 
streams in preferred stocks; it is atypical of common stocks. To soften the assumption that Xq = 0, 
one can add a constant mean to each centered self-similar X. We have investigated this possibility; 
it appears that this will require additional non-trivial tools not contained herein. Given the length 
of the current article, we have opted to leave this improvement for another work. One may, however, 
include a non-zero mean for each Xt which is proportional to this is the framework used herein 
throughout. 


1.3 Summary of main results, proof techniques, and numerics 

In [44], we studied Gaussian stochastic volatility models. The asset price process S in such a model 
satishes the following linear stochastic differential equation: 

dSt = rStdt + \Xt\StdWt, (1) 

where A is a continuous adapted Gaussian process on a filtered complete probability space (12, A, {Xt}, P), 
IF is a standard Brownian motion on (12, A, P) with respect to the filtration {Xt}, Sq = sq > 0 a.s., 
and r > 0 is the risk-free interest rate. We will assume throughout the paper that the processes 
X and W are independent. In the model in Q, the volatility is described by the absolute value of 
a continuous Gaussian process. An important special example of a Gaussian stochastic volatility 
model is the Stein-Stein model introduced in |62j , where X in 0 is an Ornstein-Uhlenbeck process. 

If Sq = Sq, the call option on S with maturity T and strike price K has price C (T, AT); this price 
equals a price Cbs {T, K-,a) in the Black-Scholes model with the volatility a depending on T and 
K. That value of a is called the implied volatility (IV) and is denoted by / (T, K). In the present 
paper, we concentrate on the behavior of C and I for small T when K is hxed; consequently, we 
typically drop the dependence of C and I on K. 

Of particular importance is the density pt of the integrated variance Yt '■= fg Xtdt. The cen¬ 
tered version of this Yt is a random variable in the second chaos of a Wiener space independent of 
W. The covariance function of X acts as a compact self-adjoint linear operator on ([0,T]), with 
non-zero eigenvalues (A^ : n = 1, 2,...) arranged in non-increasing order with repeats for multiplic¬ 
ities. This, and the corresponding eigenfunctions, are the basis for the so-called Karhunen-Loeve 
(KL) decomposition of X (see, e.g., jHESj), and of a corresponding one for Y. In any case, the 
asymptotic behavior of pr near -|-oo, which was established in |44j . depends on specific KL statis¬ 
tics, including the top eigenvalue Ai, its multiplicity ni, and the rescaled ([0, T])-orthogonal 
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projection S of Ai’s eigenspace on the mean function of X (see Theorem below). When applied 
to the case of i7-self-similar X, via the simple scaling formula priy) = , the 

behavior of pr (•) at a; —>■ +(X) translates into an expansion around T —>■ O’*" of the density pr (x) 
of the rescaled square-rooted version of Vt which is precise up to a factor (l -|- {T^)) for any 

fixed a: > 0: see asymptotic formula (181 in Theorem]^ 


Remark 1. There exist explicit formulas for the Karhunen-Loeve characteristics of various Gaussian 
processes. For Brownian motion, Brownian bridge, and OU processes, such formulas can be found 
in m- For OU bridges, one can consult [m [E], and for the Gaussian process introduced in 
[18j . the Karhunen-Loeve decomposition can be found in the same paper. Unfortunately, even for 
classical fractional Gaussian processes, e.g., fBm or fOU, the Karhunen-Loeve characteristics are 
not known. In m (see also m), Gorlay developed a powerful numerical method to approximate 
Karhunen-Loeve eigenvalues and eigenfunctions. Gorlay uses the Nystrom method associated with 
the trapezoidal integration rule combined with the Richardson-Romberg extrapolation in his work. 


The independence of W and X imply that the density Dt of St is given by a mixing formula 
(|^ involving pi via the self-similar scaling property priv) = T~^pi{T~^y). A delicate use of 
Laplace’s method then allows to translate Theorem]^ into small-T asymptotics for Dt {x) for any x 
which is “out of the money” in the context of call pricing, in the sense that the big O term depends 
on a parameter £ > 0 to allow for a: > Sq + ^ (future stock price parameter x, which stands in for 
strike price K when one computes an IV, exceeds initial stock price sq by a margin e). We find 
(Theorem that for all for x > sq -f £ 


Dt{x) = 




A,(i)--^ n 


o-iWt. 
Z 2 i 


(= 4 “) 


k^2 


Ai(l) 


Ai(I) - p/c(I) 


X X 2 log — 
So 


■»l(l)-2 


T- 


(2ff+l)ni(l) ( X 

So 


2U^l(l)r"+2 


(I + o (r^^+i)) (I + Oe ( (log 


So 


( 2 ) 


where the repeated notation (I) refers to KL elements for T = I, and where rik is the multiplicity 
of the fcth largest KL eigenvalue pk- The symbol O depends only on the covariance of A, but not 
on a; or £. The symbol depends on the covariance of X and on £, but not x. We prove formula 
([^ under the assumptions that r = 0 and the volatility process X^^'> is centered. The case where 
r > 0 and the process X^^'> is noncentered is more complicated and will be addressed in future 
publications. 

Being able to establish the precise x-behavior of the error terms above is crucial to tranposing 
the behavior of Dt (x) to the functions C and I. Specifically, we obtain the following for the 
out-of-the-money call as T —>■ 0"^ (Theorem]^ ; for AT > sqi 


C{T) = MT 


{2ff+l)(4-ni(l)) 

4 


\k) 



(3) 


where the big O above does not depend on K if it is away from sq, and the constant M is explicit 
and proportional to the constant on the right-hand side of line ([^. A nearly indentical result 
is obtained for out-of-the-money put prices P{T,K) (for 0 < K < sq) using symmetries of the 
problem (Theorem |^. 
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Ultimately, relying on a general result of Gao and Lee [35j for computing the small-time asymp¬ 
totics of IV based on those of C, we obtain in Theoremsandthat for 0 < K ^ sq, 


I{T) = 


Ai(l)-‘ 


log — 

o So 


72 


T- 


O 


__ 6H + 1 , 1 

T log - 


(4) 


where the big O is again uniform over K in any compact interval away from 0 and sq- The dominant 
factor in the expression ([^ for C, and its analogue for P, is the exponential one. In the expression 
Q for I, there is only one candidate for a dominant term. Consequently, one gets a way to estimate 
H using call or put prices or IVs away from the money as empirical statistics: 


H 


lim 

T-s-O 


log log 


C{T,K) 


lim 

T-s-O 


2 lim 
T^O 


log^ 

log log ^ 


P{T,K) 


logT 

I{T,K) 

logT 


1 

2 


1 

2 

1 

2 


where the first line holds for K > sq, the second for K < Sq, and the third holds for all K ^ sq 
(Corollaries and[^) These expressions for H do not depend on any of the model parameters 

and statistics, and are in this sense model free within the class of self-similar models. However, in 
practice, since the regime T —> 0 is limited by the ability to trade options in a liquid way sufficiently 
close to maturity, the full asymptotics in © and will typically be needed to help control the 
estimation error. 

We notice that the above asymptotics for C and I formally lose information when K = sq, since 
the expression |log(iL/so)| is zero and thus kills the dominant terms. Hence the estimators for H 
above are not longer valid in that case. We investigate this at-the-money situation in some detail. 
The delicate calculations are largely performed “by hand”. The resulting asymptotics seem to rely 
on model statistics which cannot be related to the KL elements in any simple fashion, since they 
require computing the moments ^ 1/2 and of order 1/2 and 3/2 of the non-explicit integrated 
variance’s law. As T —)■ 0, we get in Corollary that 


C(T,so) 


50/^l/2 yg+i 

7^ 



and in Theorem [TO] that 


So) — /ii/2T^ -b 


(^^ 1 / 2 ) f^3/2 ^3g_|_i ^ /rp5H+2\ 

24 \ ) ■ 


(5) 


Again, simple iL-estimators can result, which do not rely on the moments ^ 1/2 and /i 3 / 2 , such as 
Theorem ttH: 

T^O log T 

To illustrate the usage of our various asymptotic formulas numerically, we provide simulated 
stock prices, with corresponding call prices and IVs, from the self-similar volatility model, using 
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a classical a Monte-Carlo method. Using market-realistic parameter choices, we show how close 
prices and IVs are to our asymptotic formulas, noting that the fit is good in the call price case, 
and is excellent in the IV case, for time-to-maturity as large as 2 weeks. It is then not surprising 
when we show that our IV-based model-free calibration formulas for H are accurate to 2 decimal 
points up to 7 days in most cases, and 14 days in some cases. Being able to use the longest-possible 
time to maturity is important in practice because of liquidity considerations. This is all explained 
in Section [9l 

The remainder of the article is structured as follows. Some mathematical background on Gaus¬ 
sian volatility models, taken largely from [33], is in Section]^ Scaling consequences of self-similarity 
for the density of the integrated variance are provided in Section Section [^ contains the main 
asymptotic analysis of St’s density. Consequences for call, put, and IV asymptotics away from the 
money are in Sections and respectively. Sections and contain call and IV asymptotics at 
the money. The numerics in Section [^ finish this paper. 


2 Mathematical background on Gaussian stochastic volatil¬ 
ity models 

In the present section, we consider the Gaussian stochastic volatility model defined by 0. Let us 
hx the time horizon T > 0, and denote by m and K the mean function and the covariance function 
of the process X given by m{t) = E[Vt], t S [0,T] and 

K{t, s) = E [{Xt - m{t)) {Xs - m(s))], t G [0, T]^, 


respectively. It will be assumed that K{s,s) > 0 if 0 < s < T. 

The following formula is valid for the distribution density Dt of the asset price St in the Gaussian 
model described by 0: 


Dt{x) 


Vso^ - 



3 

2 


r -1 J 

f r 


i L 


log" 


X 


2 t |/2 



Pt{y)dy. 


In 0, pt is the distribution density of the random variable 


Yt = 



( 6 ) 


(7) 


The function pt is called the mixing density (see HD). The proof of formula 0 can be found in 

Ha ED. 

Applying the Karhunen-Loeve theorem to the Gaussian process .j,j, we obtain 


At = ^ ^/Ke^{t)Z^. (8) 

n=l 
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In ([^, {e„ = e„^r} is an orthonormal system of eigenfunctions of the covariance operator 

nm)= j f{s)K{t,s)ds, /eL2[0,T], 0<t<r, 

Jo,T 

and {A„ = A„(r)}, n > 1, are the corresponding eigenvalues (counting the multiplicities). The 
symbols Zn = Zn^r, n > 1, in (|^ stand for a system of iid A/’(0,1) random variables. We will 
always assume that the orthonormal system {e„} is rearranged so that 

Al — A 2 — * * * — Ayjj^ > ‘ ‘ ^ . . . 

For the sake of shortness, we introduce the following notation: 


Pi — Al, p2 — A„j+i, p3 — A„j+„2+1i ■ ■ ■ I 


5n = Sn{T) = / m{t)en{t)dt, n>l, 

Jo 

j-T 1 ni 

S = S{T) = / m{tfdt, 5 = 5[T) = 

Jo U=1 

The mixing density px is related to the density px of the integrated variance 

Yx= [ X^dt 

Jo 

as follows: 

Priy) =2Typx {Ty^) . (9) 

The next theorem, characterizing the asymptotic behavior of the density px, was established in 

mi- 

Theorem 1. If S > 0^ then the following asymptotic formula holds: 


'^1-3 


Px{x)=Cx 4 exp i/—Va: exp — 


X + O 2^^ 


as X ^ oo, where 


C = 




exp 


- 'T'^ (5^ - 'T"'" (5^ 
2Ai 


The constant A in (11) is given by 


^=n 




Al 


Al — Pk 


1^ 1 


niH-hnfc 


exp 




2 1=2 " « 


i: « 

, n=niH-hnfe-i+1 


( 10 ) 


( 11 ) 







On the other hand, for a centered Gaussian process X, we have 


Pt{x) = Cx 2 exp<j- — 


(l + 0(,T-J)) 


as X ^ oo, where 


C = 


n 


Ai 


fe =2 

The next assertion follows form Theorem [T] 
Corollary 1. The following are true: 


( 12 ) 


(13) 


1. If ni = 1, then 


Pt(x) =Cx 2 exp^6xp 
X ^1 + O 


X 

'Wi 


as X ^ oo, where C is given by (11 ). 

2. Suppose X is a centered Gaussian process with ni = 1. Then 


(14) 


Pt{x)=Cx 2 exp|-^| (^1 + 0 2^) (15) 

as X oo. 

It was established in [33] that Gaussian stochastic volatility models are risk-neutral. 

Lemma 1. In the Gaussian stochastic volatility model, the discounted asset priee process 1 1 —>■ e~'~*St 
is a {J-t}-martingale. 


3 Fractional Gaussian stochastic volatility models 

The paper |33| is mostly devoted to the extreme strike asymptotics of option pricing functions 
and the implied volatility in Gaussian stochastic volatility models. The present paper deals with 
Gaussian models, in which the volatility process is self-similar, and also with small-time asymptotic 
behavior of option pricing functions and the implied volatility in such models. 

Definition 1. Let 0 < H < 1. A stochastic process X^^'> is called i/-self-similar if for every a > 0, 
X^f^ = a^xj:^\ Here = means the equality of all finite-dimensional distributions. 

It is easy to see that if the process X^^^ is 7J-self-similar, then Xq^^ = 0. It will always 
be assumed in the sequel that the self-similar process X^^'> is stochastically continuous. For a 
Gaussian process X, the i7-self-similarity condition is expressed in terms of the covariance function 
C as follows: 

C{at, as) = a^^Cft, s), (t, s) e [0, T]^. 
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We refer the interested reader to [SlIM] for more information on self-similar stochastic processes. 
Let us consider the following Gaussian stochastic volatility model: 


dSt = rStdt + \X^t^^\StdWu So = so, 


(16) 


where sq > 0 is the initial condition for the asset price process S', W is a standard Brownian motion, 
and is a continuous iL-self-similar adapted Gaussian process. The process S characterizes the 
dynamics of the asset price in the stochastic volatility model, where the volatility is desribed by 
the absolute value of a self-similar Gaussian process. It will be assumed throughout the paper that 
the model in (16) is uncorrelated, which means that the processes X^^'> and W are independent. 


We will often suppress the parameter H in various symbols used in the paper. A popular example 
of a self-similar Gaussian process is fractional Brownian motion (see, e.g., [M])- Note that 

fractional Brownian motion is the only process that is non-trivial, self-similar, Gaussian, and has 
stationary increments. 

Exactly as in Section we will denote by pt the denstiy of the integrated variance. 


Y, = 


and by pt the density of the random variable 




ds. 


Yf. = 




ds 


(the mixing density). Since the process X^^'> is self-similar, we have Yat = Moreover, the 

following equality holds: P (Yt > y) = P (Yi > , and hence. 




(17) 


The next assertion characterizes the small-time asymptotics of the mixing density. 

Theorem 2. (i) For every a; > 0, the following asymptotic formula holds for the mixing density 


Pt in the model discribed by (16): 


Pt{x) = 2CT- 2 


exp • 


I 5(1) X 


exp< - 


2T2^^Ai(l) 


X (1 + 0, (T^)) 


(18) 


as T ^ 0, where 


C = 


A 


'ni(l) 




2V^ 


Ai(l)-i ^ 5„(1)^ 


n—l 


X exp 


Kl)-Er=l^n(l)^-E:Ll<^n(l)^ 

2Ai(l) 


(19) 
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and the constant A in {19} is given by 


^=n 


k^2 


Ai(l) 


Ai(l)-/Ofe(l) 






ni(l)H-hrafc(l) 


2 ^ Ai(l) -flfc(l) , 


(ii) If the process is centered, then 

(1 + 0, (r^)) 

as r —>■ 0, where 


C = 




A-^ n ^ 


fc=2 


Ai(l) - Pfe(l) 


fc(i) 


( 20 ) 

( 21 ) 


(Hi) If the process X^^'> is centered and ni(l) = 1, then 


Pt{x) = 2CT ^ exp i — 


2T2^^Ai(l) 


(1 + 0, (T^)) 


os r —)■ 0, where the constant C is given by (21) with ni(l) = 1. 
Proof. It follows from and 0 that 

PTix) = 2T-^^xpi (T-^^x^) . 


Since X^^) is a Gaussian process, we can use formula (|l0[). This gives 


ni (l)-3 

Pi(a:) = Cx '1 exp ■ 


/ilexpi- "" 


Ai(l) 


2Ai(l) 


X (^1 + O 2^^ 


( 22 ) 


(23) 


(24) 


(12) and (131 imply that 


as a; —)■ oo, where the constant C is given by (19). If the process X^^) is centered, then formulas 


, . „ "1(1)-^ I X 

p^{x)=Cx ^ e-PE 2 Ar(i) 


(i + o(.-i)) 


(25) 


as a; —>■ oo, where the constant C is given by (211. Now, Theorem [2| can be derived from from (23) 
(241), and ([2^. 
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4 Small-time asymptotics of the asset price density in self¬ 
similar Gaussian stochastic volatility models with cen¬ 
tered volatility. 


In this section, we restrict ourselves to the case where the process is an adapted continuous 
i7-self-similar centered Gaussian process. Recall that we assume r = 0. 

Of our interest in the present paper are asymptotic estimates of the density Dt{x) as T —>■ 0, 
which are uniform with respect to the values of a; > 0 separated from sq (away-from-the-money 
regime). Here we distinguish among two special cases. In the first case, we fix e > 0, and consider 
asymptotic expansions as t — >■ 0, which are unform with respect to a: > Sq + £■ The notation 
Oe{(j>{t,x)) as f > 0, where (j) is a. positive function of two variables, means that the O-large 
estimate holds as t > 0 uniformly with respect to a: > sq + ^- In the second case, we fix e with 
0 < e < So, and assume that 0 < a; < Sq — e. The same notation Og{4>{t,x)) will be used in the 
second case. 

Since priv) = T~^pi{T~^y), formula ^ implies that 


Dt{x) = ^r-^-2a;-5 



r 

X y exp-^ 



. 


3 

2 

fCO 

r 

X u exp ( 


-^0 

1 L 


log'f Ty^ 

2Tj/2 8 


■pi{T y)dy 


}piin)du. 


(26) 


The next assertion is one of the main results of the present paper. It characterizes the small¬ 
time asymptotic behavior of the asset price density in a Gaussian model with a centered self-similar 
volatility process. 


Theorem 3. Fix e > 0 and let x > sq + e. Then as T —>■ 0, the following asymptotic formula holds 
for the asset price density Dt in the model described by (16): 


Dt{x) = 




^ A.(l) 




2 J 


n 


X ( log - 
So 


■»l(l)-2 


(2H + l)ni(l) / X 

So 


^l(l) “ Pfc(l) 

V4 + Ai(1)t2H + 1 
2y/Xl{l)T^^2 


_ 3 

X 2 


X (1 + O (T2^+1)) ( 1 + Os ( (log |- 


(27) 


Proof. Fix a; > 0, and denote 


/•OO 

r r 

/ u ~^exp < 


Jo 

1 


log^ - 

° Sq 

2T‘^H+ly2 


+ 


T^H+1^2 


• pi{u)du 


(28) 
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It is clear from (26) that the small-time asymptotic behavior of the density Dt{x) is determined 
by that of the integral Jx{T). 

The next lemma will allow us to use Theorem to estimate the integral in 


Lemma 2. Fix a € K, & > 0, and e > 0. Let x > sq + e, and suppose f is an integrable function 
on [0,5]. Then 


[ u“exp< 


Jo 

[ 


log^f T^H+1^2 


2T-^H+ly2 
2 X 'I 


= 0, exp 


2b'^T‘^H+i 


\f{u)\dv 


as t —>■ 0. 

Proof. The lemma is trivial if a > 0. For a < 0, we have 


[ u“exp< 


Jo 

[ 


log" f- T^H+1^2 


2T-^H+ly2 


< u°‘exp 


log 


2 a; 


2T2H+ly^2 


\fiu)\du. 


\f{u)\dv 


(29) 


The following equality holds for every x4 > 0: 


A 




au 


exp<; - — 


A 

u- 


It follows that for 2x4 > —ab^, the function 


x4 

37exp<i --y 


is increasing on the interval (0, 5]. Set 


log" ^ 

x4= 


2J^2H+1 ■ 


Using (29), we obtain 


[ u“exp< 


Jo 

[ 


log" 


2T-2H+ly2 


-2 X 


log 

< 5“ exp<^~ , 


2 }j 2 rp 2 H+l 


\f{u)\du, 


\f{u)\d'u 


(30) 


provided that log" ^ > 5"T"'^+^. It is clear that the previous inequality holds for small enough 
values of T provided that cc > sq -I- £. 

Finally, Lemma follows from (30). 
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Using (23) and (25), we obtain 


Pi{y) = Ay 




exp 


2Ai(l) 




as j/ —)■ 00 , where 


)i- 


n(i) 


A = 




n '' 


fc=2 


Ai(l) - Pfc(l) 


(31) 


(32) 


It is not hard to see that Lemmaj^allows us to replace the function pi (u) in (281 by its approximation 
from (31). This gives the following: 


nOO 

J^{T)=A 

Jo 


exp < - 


log^ ^ ^ ^T^h+1 


2T‘iH+iy2 \ 


2Ai(l) 


(l + O (u du 


o.LJ^A 


2J'2H+1 


(33) 


as T —)• 0. 

To study the asymptotics of the function t i—> Jx{T) defined by (33), we consider the following 
two integrals: 


pOO 

Jx{T)=A 

Jo 


exp ■ 


log ^ /T2^+1 


2t2H+i^ 2 y 8 2Ai(1) 


du 


and 


^ poo 

XiT)=A 

Jo 


exp ■ 


log /t 2 H+i 1 


2T‘2H+iy2 ' \ 8 2Ai(l) 


du. 


Set 


Pt = 




rp2H+l 2 ^ 


2r2^+i’ "1^ 8 ' 2Ai(l)' 

Note that Pt depends on x, while 7 t does not. Then we have 


and 


= A 

[ M”i(i)-2exp(- 

IJt 2 

^ +JTU 


Jo 1 


= A 

[ u"ilil-3exp|- 

Pt 2 

^ + 7TM 


'o 1 
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du 


du. 


(34) 


(35) 


































Next, making a substitution 


.= ,|, 


we transform the previous integrals as follows: 

-i(i)-i 


and 


MT) = A[^ 
It 


J^{T)=A[^ 

It 


2 gjj-p j — -y/ 


“2 ^ 


yni(l) 3gjj.p I 


Let us denote 






dv 


dv. 


Ai(l)r2ff+i 


Then we have 


Therefore, 


\/ /3t7t = ^(r) 


log- 

So 


J^{T)=A[^ 

It 




^ u"i(i)-2exp|-z(T) 


log 


So 


dv 


and 


J,(T) =A1( ^ 
7t 


■»l(l)-2 


^ u”i(i)-3exp|-z(T) 


1 ^ 

log — 

So 


-^+v^ 


dv. 


(36) 

(37) 


(38) 


(39) 


It follows from (36) that z{T) —)■ 00 as T —)■ 0. Our next goal is to apply Laplace’s method to 
study the asymptotic behavior of the functions T 1 —>■ JxiT) and T 1 —>■ JxiT) as T —>■ 0. Note that the 
unique critical point of the function = v~‘^ + r;^ is at u = 1. Moreover, we have t/’"(l) = 8 > 0. 

We will first reduce the integrals in (38) and (39) to the integrals over the interval [0,2] and 
give an error estimate. This next assertion will be helpful. 

Lemma 3. Suppose a S K and 0 < e < Sq. Then 

J v^ exp |-\//3t7t | dv = (^exp |-3\//3 t7t|) 


os t —>■ 0. 
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Proof. Fix a small number r > 0. Then for 0 < T < Tq, we have 


V exp < — 


< c, 


■ \//?t7t ;^ + I cit'< J exp I - \/ PtItv^ | dv 

poo 

J exp (^y/P tJt - dv 

(\//3r7T - rj ^ exp |-4 (^^/;gT7T - | ■ 


The proof of Lemma is thus completed. 


Now, we are ready to apply Laplace’s method to the integrals in (381 and (39). The dependence 
of the parameter x in ( |38| ) and ( |39[ ) is very simple. This allows us to obtain uniform error estimates. 
By taking into account Lemmal^ we see that for every e > 0 and all a; > sq + e, 




7(T)=^(|) ‘ (.-(T) 


log — 
So 


exp < —2z{T) 


log- 

So 


and 


1 + Og 


z{T) 




+ Oe ( exp < —3z{T) 


1 ^ 
log — 

So 




7(T)=^(|) ‘ (.-(T) 


log — 
So 


exp < —2z{T) 


log- 

So 


(40) 


1 + Og 


z{T) 


log 


+ Oe ( exp I —3z{T) 


log- 

So 


(41) 


as T —)■ 0. Recall that the Og estimates in (40) and (41) are uniform with respect to a: > sq + e. 
Since 


J.iT) = MT) + o, (j.(r)) + Oe ( exp . 


2'J'2H+1 


as T —)■ 0, formulas (40) and (41) imply that 


= ^(|) ‘ 


log — 
So 


exp < —2z{T) 


log — 
So 


1 + 1* 

It 


Oe exp < - 


1 + Oe 


2'J'‘2H+1 


AT) 




+ Oe ( exp -3z(T) 


log — 
So 


as T —?► 0. Since for T < 1, 




( 42 ) 
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we have 

+0, (^exp|-3z(T) log ^ 

= 0, (^exp|-3z(r) log|^ I) 

= 0, (^exp|-^Ai(l)-5T-^-5 log^ 
as T —)■ 0, and therefore, 



as T —)• 0. Moreover, for all T < 1 and a; > sq + £, 



and hence 



^ Tll(l)-1 _ 1 

+ Oe(exp|-5Ai(l)-5r-^-i log^ I) 


as T —>• 0. 
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Recall that we assumed r = 0. It follows from (261 and (28) that 




2 V 2 

X exp < —2z{T) 


It 


-i(i)-i 


z{T) 


1 ^ 

log — 

So 


1 + Oe 


1 ^ 
log — 

So 


X 


log — 
So 


+ 0, (exp<j--Ai(l)-5T-"-5 


log 


So 


( 43 ) 


as T —0. 


Our next goal is to remove the last Og-term from formula (43). Analyzing the expressions in 


(43), we see that in order to prove the statement formulated above, it suffices to show that there 


exists a constant c > 0 independent of T < Tq and x > sq + e and such that 

-|Ai(l)-iT 


<cT^ 2 ( log 


So 


(2ff + l)(Tii(l)-l) „2ff + l 

8 T 4 


X ( log - 
So 


\ —2z{T) / 

X \ rr.^M+1 /, X 

— r 4 log — 

So/ V So 


The previous inequality is equivalent to the following: 


-fAi(l)-2T- 


_ (2ff + l)(Tii(l)-l) 3 / X 

< cT 8 X 2 — 

Vso 

X (log - ' 

V So 


-2RT) 


Since (42) holds, the inequality in (45) follows from the inequality 


-iAi(l)-2T-"-2 


< cT 


(2if + l)(,., (1)-1) 


X I — log — 

.■So/ V So, 


(44) 


(45) 


(46) 


Cr^e such that 


To prove the inequality in (46), we observe that for every small enough r > 0 there exists a constant 

1 

A 


C.., I - ) < log - 


. So / V So, 

for all X > So + £• Moreover, there exists Tr^e > 0 such that 




< 


So + £ 
So 




_ {2H + l)(,>l(l)-l) 

< T -8- 
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for all T < Tr^e- Now, it is clear that (46) follows from the estimate 


-H-k 


-Ai(l)-5-r)r 


> 2 + 2Ai(l)-^T-"-5-y/Ai(l)T2«+' +4 + T, 


( 47 ) 


for all T < Tr- It is not hard to see that there exist numbers t and T-r, for which the inequality in 
(47) holds. This establishes (441, and it follows that 


2 V 2 \1tJ 




z(T) log 


So 


So 


xexp -2z(T)log- 1 + oJt^ log- 


So 


(48) 


as T —>■ 0, where A is given by (32). Formula (48) will help us to characterize the asymptotic 
behavior of the function T 1 —)■ Dt{x). 

Let us assume that cc > sq + £• Then we have 


7T / V So 




(2H + 1)(t.i (1)-1) , . , 1*1(1)-! 

T 4 (1 + h) 3 


where h = -. Therefore, 


Pt 

It 


r*l(l)-l 


as T —)• 0. Moreover, 


= Ai(l) 4 log — 

V So 

(l + 0(T2-f^+i)) 

Ai(i)r2ff+i ^4 


T 


(2if + l)(ni(l)-l) 


r(T)-5 = 2 


Ai(i)r2^+i 
= i/2Ai(l)4r^ (1 + 0 (T^-f^+i)) 


and 


exp -2z(r)log- = - 


So 


So 


V4+Ai(l)T^ll + l 


(49) 


(50) 


(51) 


asT —>• 0. Next, combining (32), (48), ( |49[ ), ( |50[ ), and (51), and simplifying the resulting expressions, 
we obtain formula (27). 

This completes the proof of Theorem 
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5 Asymptotic behavor of out-of-the-money call and put pric¬ 
ing functions 

Let S be the asset price process in the model considered in ( [I^ . Define the call and the put pricing 
functions by 

C{T,K) =E[St- K]-^ and P{T,K) =E[K - St]'^ 

where T is the maturity and K is the strike price. Recall that for a Gaussian stochastic volatility 
model with r = 0, the asset price process S' is a martingale (see Lemma [^. Therefore, the put/call 
parity formula C{T, K) = P{T, K) + sq — K holds. 

In the present section, we consider the functions C and P as functions of the maturity for a 
fixed strike price, and we suppress the strike price in the symbols. Our goal is to characterize the 
asymptotic behavior as T —^ 0 of the function T C{T) for K > sq (out-of-the money call) and 
of the function T P{T) for 0 < K < sq (out-of-the-money put). 

We will first consider the call pricing function T i—>■ C'(T) with K > sq. It is known that 


C{T)= / {x — K)DT{x)dx. 


(52) 


IK 


Therefore, we can use the uniform estimate in formula (271 to characterize the small-time behavior 


of the call pricing function. Let us consider the following integrals: 


and 


h{T)= {x-K)x-i[\og- 


IK 


So 


■»l(l)-2 


exp —2z{T) log I dx 


= Sn 


Ik 


log- 
So 


r>l(l)-2 


exp S — X + 2z(T) log — > dx 


So 


-s^-^K 


IK 


log 


So 




1 “ ( d + ^^('r) I log — ^ dx 


So 


(53) 


hiT) = (x- K)x 2 ( log 


IK 


So 


■»l(l)-3 


exp —2z{T) log — \ dx 


Ik 


-s,-^K 


log- 
So 


r>l(l)-3 


exp i - 71 + 22 :(T) log — ^ dx 


So 


Ik 


1 ^ 
log — 

So 


•>l(l)-3 


exp 1 - o + 2 z(T) log — ^ dx, 


So 


(54) 


where we use the notation in (36) for the sake of shortness. 

We will next make a substitution u = (2z{T) — ^)log ^ in the integral on the second line in 
(|5^. The resulting expression is as follows: 


2ziT) - - 


^i(i) 


'(2z(r)-i)iog, 


M 3 e “du. 
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which is equal to 


(2z(r) - i 


m(1) 




0 ( -v-/ 2^ 

where the symbol F stands for the upper incomplete gamma function defined by 


POO 

r{s,x)= / v‘^~^e~'’dv. 

J X 


Making similar transformations in the other integrals in ( |53[ ) and (54), we finally obtain 

"1(1) 






So 


-s^-^k{2z{T) + \^ ' r ('^, f2z ( T ) + iVog-) 




So/ 


and 


UT) = .5 ( iziT) -1) ■ r (=i<^, ( 2.(0 - 0 log I 


-s^-^K\2z{T) + ]- 




r|!;i<^,(2.(r) + l)iog£). 


It is known that 


r(s,x) = x® ^ (1 + (s — l)x ^O {x ^)) 


as a: —)■ cx). Formula (55) can be easily derived from the recurrence relation 

r(s, a;) = (s — l)r(s — 1, a;) + x’^~^e~^ 

for the upper incomplete gamma function. It follows that 

h (T) = (^log 

^2z(r)-i +2(2z(r)-i)logf 

0 + ^ ^2(2z(r) + i)logf 

o (r^^+i^ 


2z{T) 
^K- 


= «2-(T)^-2.(T)+i 


oi(l)-2 


4z(r)2 - i 


as T —)■ 0 . Therefore, 


h{T) = (^log^) 

(i + o(r^+2)) 


oi(l)-2 


Azijf - - 


-1 


(55) 


( 56 ) 
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as T —7^ 0. Similarly, 


h{T) = (log^) 

(i + o(r^+2)) 


'^l(l)-3 


Az{Tf - 


-1 


( 57 ) 


as T —)■ 0. It is not hard to see that 


{Az{Tf =Ai(l)r2^+^ 


It follows from (561 and (57) that 


/i(r) = Ai (1)7^5 (\og- 
\ So 

(i + o(r^+5)) 


r>l(l)-2 


(I) 


2z{T) 




(58) 


as T —)■ 0. Similarly, 




/2(r) = Ai(l)7f5 (^Iog0 ' 
^i + o(r^+5)) 


(59) 


as T —0. 

The next assertion characterizes the small-time asymptotic behavior of the call pricing function. 

Theorem 4. Let K > sq. Then the following asymptotic formula holds for the call pricing function 
in the model described by (16): 


(2H + l)(4-r.i(l)) /gg 


C{T) = MT 

(14 0(T-^)) 


m 


Ai(l)“2T“"“2 


(60) 


as T ^ 0, where 


M = 






Ai(l) 


log 


K 

So 


r7l(l)-2 


n 

k^2 


Ai(l)-Pfe(l) 


(61) 


Proof. Using (27), (52), ( [53| and (54), we see that 

"1(1) 


C(T) = 


=Ai(l)- 


n 

/c=2 


Ai(i) y 

Ai(l)-pfc(l)y' 


(2H + l),il(l) 


(1 + 0 (T^^+i)) [/i(r) + o (r^hiT)^ 
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as T —7> 0. Next, (581 and (59), imply 


c(T)= 


2^T 


(=^) 


n 


Ai(l) 


Ai(l) - Pfc(l) 




('log^l ■ t ""’.;-'.,' (|)"*‘^) (n-o(r^)) 


as T —)• 0. We also have 


lAi(l)r2^^+i +4 

Ai(l)r2i/+i 


Ai(l)r2ff+ 


^ = o(r"+5) 


as T —7^ 0. Therefore, 


(|!)“'‘"=exp{-2.(r)I„g|;} 


= exp --1 


exp < 


1 /Ai(1)T2J^+i+ 4 if 


Ai(l)r2^+1 ^so 


— log — > =expi -- 


K 


jlog — 


/Ai(i)r2^+i + 4 


Ai(l)T2^+i 


2Y Ai(l)T2-ff+i ^ So 

log-1 
So 


as T —>■ 0. Using (63), we obtain 


/So\22 U) 

/so\ 

\k) 

\k) 


. 1_Ff_ 1 


Ai(1)-2T-"-2 


(i + o(r^+5)) 


( 62 ) 


(63) 


(64) 


as T —^ 0. 

Now, it is clear that Theorem follows from (62) and (64). 

The next statement allows us to recover the self-similarity index H from the asymptotics of the 
call pricing function. 


Corollary 2. Under the conditions in Theorem^ for every K > sq, 

log log CTO 1 


H = lim 
T-t-O 


log A 


(65) 


Proof. It follows from (60) that 

1 


log log i- + -"■(!)) log 1 

Ug C(T) * M 4 


/ X 1 F7 1 K _ / 2 H + l\ 

+ 2 log — + o(t—) 

So V i 


( 66 ) 
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as T —7^ 0. Hence, 


log log 


C(T) 


= log 


Ai(i)-5r 




1 ^ 
log — 

So 


+ log (^1 + 0 (^r"+5 + T^+5 log i + r"+5 0 

Ai(l)"Mog^ +0^r^+hog^^ 


= I i? + - ) log - + log 


(67) 


as T —)• 0. 

Now, it is clear that (65) follows from the previous formula. 

Next, we turn our attention to the out-of-the-money put pricing function T i— P{T) with 
0 < itr < sq. The asymptotic behavior of the put pricing function with 0 < if < sq will be 
characterized using the symmetry properties of the model in (16). In 1|41|. Lemma 9.25), several 
equivalent conditions are given for the symmetry of a stochastic volatility model. One of them is 
as follows (see (9.79) in [H]): 


= (^) Dt 


( 68 ) 


for all a; > 0 and T > 0. It is clear that for the model described by (16), the previous equality can be 
derived from formula (26). Next, using Theoremj^and (68), we establish the following proposition. 

Theorem 5. Let 0 < e < Sq and 0 < x < Sq — e. Then as T ^ 0, the following asymptotic formula 
holds for the asset price density Dt in the model described by 


Dt{x) = 






(^ 




■Ai(l)- 


7»i(i) 


n 


Ai(l) 


Ai(l)-Pfc(l) 




(2ff + l)ni(l) 


-v/4+Ai 

ry ] 2y'Vol)T’"'''2 


(?) 

2Jf+ 1 


X (1 + O (T^^+i)) (^1 + O, (log f)"’)) 


Since the model that we are studying is symmetric. 


K 


P(T,if) = -C T,^ 




K 


(69) 


(70) 


(see condition 3 in Lemma 9.25 in [41]). 

The next assertion follows from Theorem]^ and (70). 

Theorem 6. Let 0 < K < Sq. Then the following asymptotic formula holds for the put pricing 
function in the model described by {16): 


P{T) = MT 


(2ff + l)(4-r.i(l)) 


(140(T^)) 


Ai(l)-2T- 


(71) 
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as r —>■ 0, where the constant M is given by 


M = 


{soK)l 


Z 2 i 




A.(l) 


4-n-,(l) 


(log 


r^l(l)-2 


OO 


xn 


^i(i) 

Ai(l) - Pfc(l) 


2 


(72) 


Next, using the same reasoning as in 
Corollary 3. Under the conditions in 

H = 


the proof of Corollaryj^ we obtain the following statement. 
Theorem^ for every 0 < K < sq, 
log log 1 


lim 

T-fO 


log A 


(73) 


6 Asymptotic behavior of the implied volatility 


Theorems and characterize the small-time behavior of the call and put pricing functions in a 
stochastic volatility model with centered Gaussian self-similar volatility. In the present section, we 
study the small-time behavior of the implied volatility in such a model. We will use some of the 
results obtained by Gao and Lee in . Gao and Lee establish certain asymptotic relations between 
the implied volatility and the call pricing function under very general conditions. They consider 
various asymptotic regimes, e.g., the extreme strike, the small/large time, or mixed regimes. Of our 
interest is formula (7.11) in Gorollary 7.3 in |35j . providing an asymptotic formula characterizing 
the small-time asymptotic behavior of the implied volatility in terms of the call pricing function. 
It follows from this formula that UK ^ sq, then 


K) = 


log — 

o So 


log 


C{T,K) 


as T —)■ 0. Therefore, 


1 - 1-0 


log 

1 1 

iOg c(T,K) 





iOg C{T,K) 



I{T,K) = 


log — 

® So 


2T 


log 


O 


log 


log 


C(T,K) 


C{T,K) 


Vt 


log 


C(T,K) 


(74) 


as T —0. 


The following assertion can be derived from (601 and (74). 

Theorem 7. Let K > sq. Then the following asymptotic formula holds for the implied volatility 
in the model described by (16): 


I{T) = 


Ai(i)^yi^ 


V2 


T- 


+ 0 T 


log 


T 


(75) 


as T —0. 
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Proof. It follows from (66) and (67) that 


C{T) ^ 


-H-l 


and 


log log 


C{T) ~ 


as T —?► 0. Moreover, the mean value theorem implies that 


log 


1 


C(T) 


= (Ai(i)-5r-"-hog^^ Vo('r'^iog^ 


. . .L , K\ ^ 2H + 1 ( 6H + 3 1 

= Ai(l)'‘ log—) T 4 +o(t 4 log- 


so 


T 


< 


as T —)■ 0. Now it is not hard to see that (75) follows from (74) and the previous formulas. 

Remark 2. Assume K > sq . It follows from Theorem that if the Hurst index satisfies 0 < H 
then the implied volatility T I{K,T) is singular at T = 0, and it behaves near zero like the 
function T i— T 4 . For standard Brownian motion, H and we have 


2 ’ 


lim /(AT, T) = 


Ai(l)iyiogf 


T-s-o ' ' ' y /2 


Finally, for | < i7 < 1, the implied volatility T 1 —>■ I{K,T) tends to zero like the function T 


^2H-1 


The next statement is a corollary to Theoremj^ It provides a representation of the self-similarity 
index in terms of the implied volatility. 


Corollary 4. Let K > Sq. Then the following equality holds: 

log 1 


H = 2 lim 


I(T,K) 


(76) 


T^o log ^ 2 

In the case where 0 < K < sq, Theorem 0 Corollary m and the symmetry condition 

(see mi, Lemma 9.25) imply the following assertions. 

Theorem 8. Let 0 < K < sq- Then the following asymptotic formula holds for the implied volatility 
in the model described by (16): 


, , Ai (1) 4 ^/log ^ 2H-1 ( en+i 1 

/(T) = + O ( t ^ log - 


( 77 ) 


as T —0. 


Corollary 5. Let 0 < K < sq- Then equality (76) holds for the self-similarity index H. 
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7 At-the-money options 

In this section, we consider a stochastic volatlity model, in which the volatility process is an 

adapted 77-self-similar Gaussian process. As before, we assume r = 0. Let us also suppose K = sq 
(at-the-money case). Note that here we do not assume that the volatility process is centered. 


Using (26) and the formula 


poo 

C{T,K)= {x - K)DTix)dx, 
Jk 


we obtain the following equalities for the at-the-money call: 


C'(r,so) = ^T-^-^ r 

V 27r Jo 

X / (a; — exp <— 

So ( 


u ^exp < — 


T^H+1^2 


Pi {u)du 


2T2H+ly^2 


dx 




X 2 exp ■ 


u ^exp < — 


T^H+1^2 


Pi(u)du 


log" 


2t‘^h+i^2 


dx — So X 2 exp < — 


log" 


2 T‘^H+ly 2 


dx 


It follows from the previous formula that 


u ^ exp 


T2H+1^2 


Pl{u) 


X [^i('r, u) - ^2{T,u)]du, 


where 


and 


poo 

j ?/"5exp< 

f log" y \ 

( 2r2^+iu2 j 

r i i 

log" y \ 

1 1 / = exp <1 

2T^h+i^2 ^ 



(78) 

dy 

(79) 

dy. 

(80) 


Our next goal is to estimate the functions $1 and $2 defined in (79) and (|80|. We have 

,,2 


$i(r,M) = ^ expj- 


|y2H+1^2|^oo 


= exp 
= exp<j 


w w 

2 T-^H+ly 2 ~ ^ 

f 1 


dw 


exp ■ 


r t^h+ 1^2 


_ rpH+\ 


8 


/_lT2ff + l„2 


2 T‘iH + ly 2 

exp<^ - 


1 




uexp ■ 


2 T-^H+iy 2 

2 


T^H+1^2 


dz 


dw 
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Similarly, 


$2 {T, u) = 2 u exp 


T^H+1^2 

8 



exp 



dy. 


Therefore 


<&! (T, m) - $2 (T, u) = 2r^+ 2 u exp 


T^H+1^2 


\T ^2-i 


exp< - 


y 


dy. 


(81) 


The next lemma will be useful in the sequel. It will allow us to estimate the integral in (811. 
Lemma 4. Let 0 < a < 1. Then the following inequalities are valid: 


6 

On the other hand, if a> 1, then 

1 


o —;r ^ / 6xp } dy < a —^t + tx- 


\/2 a 


exp 


a 

'T 


< / exp - 


< 


6 40 


dy 


(82) 


Jo 

a 

a^ + l 


exp 


a 

'T 


(83) 


Proof. The inequalities in (82) can be established using the Taylor expansion with two and 
three terms. 


To prove the estimates in (83), we use the following known inequalities: 


+ 1 


exp< - 


< 


exp < - 


y 


dy < — exp ■( — 

X ' 


(84) 


for all X > 0. The previous inequalities follow from stronger estimates formulated in [T], 7.1.13. 
Now, (83) can be derived from (84) and the equality 


exp <i Hy = ^ - 


exp<{ -y \dy. 


This completes the proof of Lemma 

The next assertion provides estimates for the at-the-money call. 
Theorem 9. The following inequalities are true for every T > 0; 

Ui{T)<C(T,so)<U2{T), 


where 


nOO nQ 

U,{T) = / p,(u)udu - ^ 

v 27 r Jo 24 v 27 r Jo 


25n 




Pi{u)u^du 




dv 
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and 


U2{T) = 

V 27r 


Pi{u)udu - 

24v^ 


+ 


So 


Pi{u)u^du 
V 


6A0V^ 




V 

7^ ~ 2 48 ~ 1^ 


2 v 


- exp 


dv 


Proof. It follows from (78), (81) and Lemma|^that 

2 

'^O / 7 


C{T,so) < 


727r . 
1 


T^+hu- — 
24 


Pi{u) 

1 

(340 




du 


2 so 

7^, 

So 

' 

2 so 

7^ 

2 so 

7^, 


Pl{u) 


^H+i 

poo 


^/7T 

2T^+^u f 

T^H+1^2Y 

[72 

T^h+ 1^2 + 4 t 



du 


Pi{u) 


/: 


1—1 y5//+|^5 


T^+5m-+ 

24 640 


du 


r"+5 

/•OO 


Pl{u) 

^T^+ku- 
2 48 

Pl{u) 

7 ^ 2T^+hu 

72 T^H+iy 2 4 


3 Q 
2 li'^ 


exp 


1280 




dw 


y 2 ff+ly 2 


du. 


(85) 


and 


C(T,so) > 


So 

7^ 


Pl(u) 


T^+iu- —T3^+57 
24 


2s 


0 


7^. 

So 

72 

2 so 


Pi (7 


= w.l *<“> 


72 r^+5u 

1 


exp< - 


T^H+1^2 


du 

du 


24 


du 


7^ 

2 so 

7^ 


Pi (7 
Pi (7 




48 


.72 T^+i 


exp < - 


du 


T^H+1^2 


du. 


( 86 ) 


Now, it is not hard to see, making the substitution v = T^^^u, that Theorem follows from 
(^ and 

The next statement characterizes the small-time asymptotic behavior of the at-the-money call 
pricing function in a Gaussian self-similar stochastic volatility model. 
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Corollary 6. The following formula holds as T ^ 0; 


where 


and 


C{T, So) = CiT^+5 - + Q , 

) 

Pi{u)u^ du 


Cl = 


C2 = 


So 

j 

So 

24v^. 


Pi{u)u^ du. 


(87) 

( 88 ) 

(89) 


Proof. For a centered volatility process X, we will use formula @). In the case of a noncentered 
volatility process X, we need the following formula: 


Pi (x) = 2C'x 2 exp ■ 


I S{1) 

Ai(l) 


X > exp 


2Ai(l) 


X (l + 0 (x-i)) 


(90) 


as X —>■ oo, where the constant C is given by Formula ( [90| now derives easily from ( [2^ and 

It follows from Theorem [9] that 

C{T,so)-UiiT) < U 2 iT)-UiiT) 

So 


< 


POO 

T^^+l / p^{u)u^du 

Jo 


640v^ 

2 so 




-exp<^ - — 

V 


2v 


exp i —JT 


1280 


dv. 


(91) 


Let us next suppose the process X is centered. Then, using (311, we see that for x > 2 and for 
sufficiently small values of T, 


1 




< a 


< a 






1 


1 




exp 


„,2 


exp 


2Ai(l)T2^+i 


< a- 


1 


< a exp < — 


exp 


4Ai(l)r2^+i 

1 

'2Ai(l)r2^+i 

1 

exp 


exp 


8Ai(l) 


^ 4Ai(l)T2^+i 
Here a > 0 is a constant that may change from line to line. 


8Ai(l)J ■ 


(92) 


30 































Now assume the process X is noncentered. Then for z; > 2 and for sufficiently small T, 


1 ^ I V \ IV 

rPi I I < a 


T^+h \T^+ 


T^+h 




<5(1) 


T^+h. 1 Y Ai(l) 


exp< - 


< a 


2 Ai(l)T2^+i 


1 


T^+5 


exp < - 


4Ai(l)r2^+i 


< a exp < — 


4Ai(l)T2^^+i; 8Ai(l);' 


Finally, taking into account (p^, (921, and (p^, we obtain 


C'(T,so) - Uy{T) = O 


as T —)■ 0. Now, it is not hard to see, using the definition of ?7i, (92), and (94) that 


C{T,sq) = 6iT^+5 - + o 


where 


and 


bi = 


bo — 


So 

v^. 


So 


Pi(u)udu 


Pi{u)v?du. 


24 -\/ 27 r Jo 

Finally, using the equality pi(u) = 2upi(v?'), we obtain bi = Ci for z = 1, 2. 
This completes the proof of Corollary 


(93) 


(94) 


8 Implied volatility in at-the-money regime 


The Black-Scholes call pricing function for r = 0 and iF = sq is given by 

gy/T 


CBs{T,So,cr) = 


So 




dy = So—i= 


. crVr 
2\/2 


dx. 


/o 


Hence, 


CBsiT,so,(j) = So erf ^ 


(95) 


where erf is the error function defined by erf(u) = e ^ dx. The error function is a strictly 

increasing continuous function from [0,cx)) onto [0,1). Its inverse function is denoted by erf~^. It 
is known that the inverse error function has the following Maclorin’s expansion: 


7tt^ 


err\z) = y- z + -z^ + —z 


12 


480 


0 < z < I 


(96) 
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(see []). It follows from the definition of the implied volatility that 

CBsiT,so,I{T,so)) = C{T,so). 


Therefore, (95) implies 


Next, using (96), we obtain 


IiT,so) 


Vt \ So 


I{T,so) 


r C{T, So) 
Vt [ So 


TT C{T,sof 
12 sg 


+ 0(C(T,soV) 


(97) 


as T —)• 0. 

Now, we are ready to characterize the small-time asymptotic behavior of the implied volatility 
in at-the-money regime. 

Theorem 10. The following asymptotic formula holds as T ^ 0; 

pOO 

I{T,so)=T^ / pi{u)u^du 

Jo 


I rp3H-\-l ^ 

24 

-f O (T^^+^) . 


l‘Oo ^ \ 3 poo ^ 

/ pi{u)u^du\ — / pi(u)u^du 


(98) 


Proof. Our first goal is to obtain an asymp toti c fo rmu la for the implied volatility with error 
term of the order O by using formula (87) in (97). Following this plan, we obtain 

/(T, so) = ^ (ciT^+s - C2T3^+i + O 

■ (^c^T^+i-C2T^^+i +o(T^^+i)y+ 0{T^^+V 


12s?,Vt 


'VJVci rpH I f ^ VTkC2 \ rpZH+1 I Q (ji5H+2\ 

So y 12si) So y ^ ^ 


(99) 


as T —0. Now, it is not difficult to see that formula (98) follows from (88), (89), and (99). 
This completes the proof of Theorem 


Remark 3. It is clear that the following formulas are valid for the integrals in (98): 

poo ^ 

Mi/2 := / pi{u)u2du = E 


Xlds 


and 


M3/2 



r 1 

/pi \ 2 

3 

/ pi{u)u"^du = E 

( / X^ds) 

Jo 

\Jo ) 
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Theorem 10 allows us to recover the self-similarity index H knowing the small-time behavior of 
the at-the-money implied volatility. 


Theorem 11. The following formula holds: 


H = lim 
T-J-O 


logi 


9 Numerical illustration 


To illustrate the numerical potential of our asymptotic formulas in practice, we finish this arti¬ 
cle with a brief section comparing exact (Monte-Carlo-simulated) option prices and IVs with the 
asymptotics we have derived. Formulas such as Q can be used to calibrate various parameters 
which might be linked explicitly or empirically to Ai (1), assuming H is known. We refer to the 
numerics in our prior work in |44j for details on what can be done, leaving to the interested reader 
any details of how to translate the ideas therein which are for extreme strike asymptotics to the 
small time case. [H] also contains a description of how to simulate the fBm-driven models of inter¬ 
est to us, for Monte-Carlo purposes, as alluded to in Remark we do not repeat this information 
here. 

Our results in the at-the-money case are presumably harder to exploit along these lines because 
they depend on moment statistics fJ. 1/2 and /i 3/2 (Remark [^, which are not explicitly related to 
model parameters. An exception to this observation is in the case of models with a volatility scale 
parameter a, by which we mean that one replaces model Q with 

dSt = rStdt + a \Xt\StdWt- ( 100 ) 


Here the parameter a is rather inoccuous since, by self-similarity of |A|, this cr can be absorbed 
as a linear time change, but it represents a convenient parameter for tuning a model to realistic 
time-scales and volatility levels. We will use this device in this section. In particular, at the money, 
it is easy to see from Theorem that one has 


where /ii /2 and 11^/2 are given in Remark Thus at-the-money IV asymptotics can be used to 
calibrate tr in model (100). We do not comment on this further herein. 

Instead, we provide a numerical analysis of our results’ use in iJ’s calibration. Indeed, the 
reference [H] contains an effort to calibrate H itself, when other parameters have been estimated 
by other means, but left some stones unturned. We found therein that H calibration can be 
relatively successful in some cases in practice, though this is not necessarily backed up by any 
asymptotic theory. In this section we show instead how model-free results such as Corollary and 
Theorem E] provide excellent calibration of H in many cases. We choose to present this in the 
at-the-money case for two reasons. First, it illustrates the model-free framework, since the results 
we obtain are not sensitive to the values of ^ 1^2 and p. 3 / 2 . Second, in practice, liquidity is low for 
options away from the money near maturity, which all but dictates the use of at-the-money IV. 

The setup we use is that of model (100) with X = fBm, r = 0, and a = 3. The choice of a 
is tailored to provide a realistic volatility level after 1 or 2 weeks, with time measured in years. 
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Specifically, a practicioner may simply select the desired magnitude of tr by matching it to the mean 
magnitude of volatility in (100) via the formula 


E[a\Xt\]=<7t^y^. 

For example, with H = 0.6 and cr = 3 we get E [a |Xt|] « 0.22 after one week {t = 7/365 « 0.019,) 
and E [a |Xt|] « 0.34 after one week {t = 14/365 « 0.038), which could represent a realistic scenario 
for a volatile short-term bond market. Values of a closer to unity result in much smaller volatility 
values near maturity; these allow for an extremely sharp fit between theoretical call and IV values 
and our asymptotics, but would typically be unrealistically small, hence our choice of cr = 3. 

a first question might be whether it would not be sufficient to use 


Before using Theorem 11 


an asymptotic theory for call prices to estimate parameters. The use of IV over option prices has 
been advocated in many articles, including many of the ones cited herein, but the question is still 
legitimate since one rarely sees evidence in the literature that this is indeed preferable in practice. 
The two images in Figure 1 compare the fit between our asymptotic formulas (Corollary and 
Theorem 101 and exact (simulated) call and IV values for times from 1 day to 2 weeks. 



time time 


Figure 1: Call (left) and IV (right) with a = 0.3, t £ [1 day; 2 weeks], H = 0.51 

We chose the extreme case H = 0.51 because, as it turns out, the asymptotics’ accuracy increase 
as H increases. We see from the above that the IV asymptotics are accurate at a roughly 5%-error 
level for more than 10 days, and remains fairly accurate up to 2 weeks, while the call asymptotics 
are only accurate at a 5%-error level for 2 days, and deteriorate significantly thereafter. Other 
values of H show similar pictures. The choice to use IV over call prices for calibration purposes in 
small time is clear. This can of course be verified rigorously on our formulas since our coefficients 
can be computed numerically as well; this is omitted from our study. The four pictures in Figure 
2 show the extremely sharp fit of IV asymptotics over two weeks as H increases, as we mentioned. 

Since liquidity decreases as time to maturity decreases, it is desirable to use the largest possible 
time to such that the relative error in IV approximation does not exceed a given error level, say 
1% which would be a high level of accuracy. The table below give an idea of what this means in 
practice, by computing to for a 1% level in the above realistic cases: with 

f simulated IV (t) — asymptotic IV (t) „ „ 1 

to = max t : - . - — < 0.01 \ 

\ simulated IV (t) J 
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Figure 2: IV with a = 0.3, t S [1 day; 2 weeks], H = 0.55,0.60 (top left and right), H = 0.75,0.85 
(bottom left and right). 


we find : 


H 

0.51 

0.55 

0.60 

0.75 

0.85 

to in days 

2.3 

4.6 

10.4 

14 

14 


These values of to could be considered as rather conservative, due to the choice of 1% accuracy; 
practitioners may decide to choose a slightly more liberal level. This is evident from the last tables 
below, in which we show the result of the calibration of H from exact (simulated) option prices, 
via Theorem HH 

T = 1 day 


H used in simulation 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 

H calibrated from IV via Theorem 

11 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 


T = 2 days 


H used in simulation 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 

H calibrated from IV via Theorem 

11 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 
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T = 7 days 


H used in simulation 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 

H calibrated from IV via Theorem 

11 

0.50 

0.50 

0.55 

0.60 

0.75 

0.85 


T = 14 days 


H used in simulation 

0.50 

0.51 

0.55 

0.60 

0.75 

0.85 

H calibrated from IV via Theorem 

11 

0.50 

0.50 

0.54 

0.59 

0.75 

0.85 


In all cases, even with a 14-day time to maturity, the error in iJ-cahbration is no greater than 
one hundredth (less than 2% relative error). The only difficulty we experience appears to be in 
differentiating between a model with Brownian scaling (iJ = 0.50, no memory in the volatility) 
and a model with H > 0.50, except for the very short times to maturity t = 1, 2 days. If liquidity 
at those levels is adequate, as it may be in heavily traded bond markets, then our calibration can 
be used with such short horizons. Otherwise a maturity of one week is preferable, particularly for 
self-similarity indices which are not too close to 0.50. A maturity of two weeks will work in all cases 
for scenarios where one is satisfied with a possible error of one hundredth on H calibration; this 
could be a realistic accuracy level for many users of stochastic volatiltiy models who are currently 
not using any self-similarity or long-memory assumptions. 
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